Simulation of a quantum phase transition of polaritons with trapped ions 
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We present a novel system for the simulation of quantum phase transitions of collective internal 
qubit and phononic states with a linear crystal of trapped ions. The laser-ion interaction creates an 
energy gap in the excitation spectrum, which induces an effective phonon-phonon repulsion and a 
Jaynes-Cummings-Hubbard interaction. This system shows features equivalent to phase transitions 
of polaritons in coupled cavity arrays. Trapped ions allow for easy tunabilty of the hopping frequency 
by adjusting the axial trapping frequency, and the phonon-phonon repulsion via the laser detuning 
and intensity. We propose an experimental protocol to access all observables of the system, which 
allows one to obtain signatures of the quantum phase transitions even with a small number of ions. 
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Trapped ions are among the most promising physical 
systems for implementing quantum computation [l| and 
quantum simulation [2J. Long coherence times and indi- 
vidual addressing allow for the experimental implemen- 
tation of quantum gates and quantum computing proto- 
cols such as the Deutsch-Josza algorithm, teleportation, 
quantum error correction, quantum Fourier transforma- 
tion and Grover search [3|. Quantum simulation could 
be performed in the future by large-scale quantum com- 
putation [4]. With the currently available technology, 
tailored Hamiltonians can be modeled with trapped ions 
to simulate mesoscopic Bose- Hubbard systems @, spin- 
boson systems [6*] and spin systems [?} . Recently, a quan- 
tum magnet consisting of two spins has been successfully 
simulated experimentally with two trapped ions Q. 

In this Letter, we propose a physical implementation 
of the Jaynes-Cummings-Hubbard (JCH) model using 
trapped ions. The JCH model was proposed in the con- 
text of an array of coupled cavities, each containing a 
single two-state atom and a photon Q . Such a system is 
described by the combination of two well-known physical 
models: the Hubbard model [13, Eij which describes the 
interaction and hopping of bosons in an optical lattice, 
and the Jaynes-Cummings model, which describes the in- 
teraction of an atom with a quantum field . The JCH 
model predicts a quantum phase transition of polaritons, 
which are collective photonic and atomic excitations. We 
shall show that the laser-driven ion chain in a linear Paul 
trap is described by a JCH Hamiltonian, wherein the ions 
and the phonons correspond, respectively, to the atoms 
and the photons in a coupled cavity array. Fig. [TJ As 
in 0], the position-dependent energy and the non-local 
hopping frequency of the phonons is controlled by the 
trapping frequencies, while the effective on-site repulsion 
is provided by the interaction of the phonons with the 
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FIG. 1: (Color online), a) Coupled cavities each contain- 
ing photons and single two-state atoms. Intercavity hopping 
is provided by an optical fibre. The strong coupling between 
the atoms and the photons leads to an effective photon-photon 
repulsion, b) All ions are simultaneously interacting with a 
traveling wave in the radial direction. The laser-ion interac- 
tion creates an effective on-site interaction between the local 
phonons. The phonon hopping appears due to the Coulomb 
interaction and can be adjusted by the mutual distance of the 
ions. 



internal states of the ions and can be adjusted by the 
parameters of an external laser field, namely the Rabi 
frequency and the detuning. This on-site interaction is 
analogous to the photon blockade (photon-photon repul- 
sion), where the strong atom-cavity coupling prevents the 
entering of additional photons into the optical cavity [l^ . 
We shall show that many-body effects appear as a quan- 
tum phase transition between a localized Mott insulator 
(MI) and delocalized superfluid state (SF) of the com- 
posite phononic and internal (qubit) states of the ions. 
Due to the collective nature of the excitations we distin- 
guish between collective qubit and phononic SF and MI 
phases, and the pure phononic SF phase, similar to the 
effects predicted in [1J| for coupled cavity arrays. 

Consider a chain of N ions confined in a linear 
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Paul trap along the z axis with trap frequencies Wg 
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{q = X, y, z), where the radial trap frequencies are much 
larger than the axial trap frequency {ujx.y ^ w^), so that 
the ions are arranged in a linear configuration and oc- 
cupy equilibrium positions along the z axis. Making 
a Taylor expansion around the equilibrium position, and 
neglecting , y^, zx^, zy^ and higher order terms, the 
Hamiltonian in the radial direction x reads 
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Here pk is the momentum operator, M is the ion mass, Xk 
is the position operator of the fcth ion about its equilib- 
rium position Uk — z^/l^ where I — {e^ /'ineoMuj'^y^'^ 
is the length scale, e is the charge of the ion, and Eq is 
the permittivity in free space. In the Hamiltonian ([1]) the 
motion in the radial direction is decoupled from the axial 
motion. In terms of the normal modes Wp, the Hamilto- 
nian ^ reads = fiJ2p=i'^p(^lo-p+^)- Here a J and cip 
are the phonon creation and annihilation operators of the 
pth collective phonon mode. However, if Xk and pk are 
written in terms of local creation aj^ and annihilation dk 
phonon operators, so that Xk = ^h/2AIuJx (a^ +afe) and 
Pk ~ iy/hMuJx/2 (a^ — ftfe), the Hamiltonian ([T|) reads 
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where we have neglected higher-order (energy non- 
conserving) terms. The phonons are trapped with a 
position-dependent frequency 
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where a = ujz/^xi and they may hop between different 
ions, with non-local hopping strengths 
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derived from the long-range Coulomb interaction [s^ . 

The collective and local creation and annihilation op- 
erators are connected by the Bogoliubov transformation. 
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which preserves the commutation relation, [dfe , dj„] — 
6km- Here 9p = — ^ ln7p, with 7^ = 1 + Q!^(1 — Ap)/2 and 
\p are the eigenvalues, with eigenvectors b^P\ of the ma- 



trix Ai. 
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Using Eq. (O one finds that the pth collective phonon 



Here \2nk) {k = 1,2...N) is the local Fock state with 
2nk phonons. For the center-of-mass phonon mode we 
have p — 1 and cosh 9p — 1; hence the collective 
ground state |0i) is a product of local ground states, 
|0i) = |0i) . . . IOat). For a sufficiently small number of 
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Since the 



lowest-energy collective vibrational mode in the radial 
direction is the highest mode p = N, we find that the 
superfluid ground state of the Hamiltonian ^ is 
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Here we have assumed the commensurate case where the 
number of ions is equal to the number of phonons. We 
find that the ratio between the average number of local 
phonons in the ground state is given by the square of the 
oscillation amplitudes: {fik) / {fim) — {^^j^^ l^m^Y , where 
fik = aj^dfc {k = 1,2,..., N) is the local phonon number 
operator. 

We shall show that the laser-ion interaction induces 
an effective repulsion between the local phonons. This 
interaction provides the phase transition from phononic 
SF state to composite SF and MI phases of the joint 
phononic and qubit excitations. Consider ion qubits with 
a transition frequency loq, which interact along the ra- 
dial direction with a common traveling-wave laser light 
addressing the whole ion chain with frequency wl- The 
Hamiltonian of the system after the optical rotating- wave 
approximation is given by 
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Here a^T = \ek) {gk\ and CTj, = \gk) {ek\ are the spin flip op- 
erators, |efe) and \gk) are the qubit states of the fcth ion, 
is the real- valued Rabi frequency, 5 = uji^ — loq is the 
laser detuning, and rj — |k|xo is the Lamb-Dicke param- 
eter, with k the laser wave vector, and xq = ^/h/2Muj^ 
is the spread of the ground-state wave function. The 
Hamiltonian, after transforming into the interaction pic- 
ture by the unitary transformation U = Q^Hat/n^ with 

ffo = -f>^xY.k=i0^k°'k + f^^Ylk=i |efe)(efcl> in the Lamb- 
Dicke limit and after the vibrational rotating-wave ap- 
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FIG. 2: (Color online) Chemical potential /i as a function 
of the laser detuning A. Each curve corresponds to ^J.{n) — 
Eg{n + 1) — Eg{n), where Eg is the ground energy of the 
single ion Hamiltonian Q for zero hopping. Here n is the 
number of excitations per ion. The Mott lobes are defined as 
Sfik = ii{k + 1) - fi{k), fc = 0, 1, 2 . . . 
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where Hi = W HU — ihW dtU ■ We assume that the laser 
is tuned near the red motional sideband S — —Ux — A, 
with a small detuning A (A ^ lOx). The coupling 
between the internal qubit and local phonon states is 
g = r/f2. The Hamiltonian ^ is valid when tkm,g ^ ^x, 
which ensures that higher terms, which violate the con- 
servation of the total number of excitations, can be ne- 
glected. The first three terms in Eq. ^ describe the 
Jaynes-Cummings model. The first two terms correspond 
to the energies of the local phonons and the ions, while 
the third term describes the laser-ion interaction. The 
fourth term in Eq. ^ describes the non-local hopping of 
phonons between different ions and allows the compari- 
son to Hubbard systems. The last term —^N is the usual 
chemical potential. 

The Hamiltonian ^ commutes with the total excita- 
tion operator N = J2k=i ^k, hence the total number of 
excitations is conserved. Here Nk = al.ak + \ek){ek\ is the 
number operator of the total qubit and phononic excita- 
tions at the kth site. The effective on-site interaction is 
provided by the interaction of phonons and qubit states 
at each site. The strength of the on-site interaction de- 
pends on the external parameters, such as the Rabi fre- 
quency n and the laser detuning A. This interaction 
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FIG. 3: (Color online) Total (qubit -|-phonon) variance 'DNk 
(top) and the qubit variances VNa,k (k = 1,2,3) (bottom) 
for a chain of five ions with five excitations as a function of 
the laser detuning A for fixed hopping t = 0.3g. Negative 
values of A correspond to blue detuning with respect to the 
red-sideband transition. 



creates an energy gap, which prevents the absorption of 
additional phonons by each ion. In Fig. [2] we plot the 
chemical potential fi, which counts the energy to add a 
extra excitation into the system, as a function of the laser 
detuning A for zero hopping. The Mott lobes Sfik > for 
sufficiently large number of phonons tend to zero, which 
show that the on-site interaction is vanishing [9] . 

We describe the quantum phase transition between the 

MI and SF states by the variance VNk = \J {Nl) - {Nk)"^ 

of the number operator Nk with respect to the ground 
state of the Hamiltonian ^ for fixed number of excita- 
tions . If the on-site interaction between the phonons 
dominates the hopping, the ground state wave function 
is a product of local qubit and phononic states for each 
site with a fixed number of excitations. Hence in the 
MI state, the variance DNk for any k vanishes. When 
the hopping term dominates the on-site interaction, then 
the ground state consists of a superposition of qubit and 
phononic states with delocalized excitations over the en- 
tire chain. In this state the variance VNk at each site 
(i.e. each ion) is non-zero. 

Figure [3] (top) shows the variance VNk {k = 1,2,3) 
for a chain of five ions with five collective excitations 
versus the laser detuning A for fixed hopping frequency 
t = aLOzl^ calculated by an exact diagonalization of the 
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Hamiltonian ([9]) . Due to the symmetry of the trap with 
respect to the center it is not necessary to plot the phase 
diagrams for ions 4 and 5. For sufficiently large negative 
detuning A, there exists an energy gap, which prevents 
the absorption of additional phonons. Hence, the system 
is in the MI phase, where the qubit and phononic excita- 
tion are localized. When the detuning A increases, the 
energy gap decreases and the system makes a phase tran- 
sition to the SF phase. The phase transition is stronger 
for the ions near the center of the trap due to stronger 
Coulomb interaction and increased hopping strengths, 
and weaker at the ends of the ion chain. The comparison 
between the variance at the different sites demonstrates 
two characteristic features. Firstly, there is a range of de- 
tunings where the qubit and phononic excitations at ion 
1 (end of the chain) are predominantly in a MI phase, 
whereas the other ions are in the SF phase. Secondly, 
there is a broad range of A along which the joint excita- 
tions at all ions are in the SF phase. 

Figure [3] (bottom) shows the variance of the qubit ex- 
citations VNaM with Na,k = \ek){ek\, [k = 1,2,3) [H. 
This allows us to distinguish the following phases: in the 
region of large negative detuning A the collective and 
the qubit variances are small, indicating that the sys- 
tem is in the qubit MI phase. Increasing the detuning, 
the collective variance stays small but the qubit variance 
increases, which shows that the system is indeed in a col- 
lective MI phase. Approaching A = the system makes 
a phase transitions into the collective qubit and phononic 
SF phase as now both collective and qubit variance arc 
large. Finally, for sufficiently large positive detuning the 
qubit variance decreases but the collective variance stays 
large, which shows that the system is in the phononic 
superfluid phase. 

The experiment is started by initializing the ion chain 
with N phonons in the lowest energy radial mode. To 
avoid off resonant excitation of unwanted radial modes, a 
could be increased temporarily. Then the coupling laser 
is switched on. The experimental proof for the phase 
transition can be carried out by local measurements re- 
peated for each ion which should be performed faster 
than the hopping time. In the case of ^°Ca+ ions the 
qubit states could be represented by the ground state S1/2 



and the metastable state D5/2- The laser, which creates 
the phonon-phonon repulsion, would be detuned to the 
red sideband of the quadrupole qubit transition between 
the two states. Then the readout could be performed by 
scattering photons on the dipole transition S1/2 ~>Pi/2- 
This would lead to momentum recoil and changes of the 
phononic excitation, but to circumvent this we have to 
perform a measurement of the qubit states, which does 
not affect the phononic state as in the following steps. 
1) Make a random guess for the qubit excitation. 2) 
If the guess was the S1/2 ground state, then swap the 
population S1/2 "^05/2 by a carrier tt pulse leaving the 
phononic excitation unchanged. 3) Now expose the ion 
to laser light on the dipole transition. 4) If the ion scat- 
ters light the guess was wrong and we have to discard 
the measurement and restart, otherwise the initial guess 
was right and we transfer the qubit excitation back to the 
Si/2 ground state by another carrier tt pulse, then drive 
Rabi oscillations on the red sideband by perpendicular 
Raman light beams with the difference momentum vec- 
tor in the radial direction. 6) The phononic population 
can now be extracted by a Fourier analysis of the Rabi 
oscillations. 

In conclusion, we have proposed a novel implementa- 
tion of the JCH model by trapped ions simulating po- 
laritonic phase transitions in coupled cavity arrays. The 
system shows a Mott insulator to superfluid phase transi- 
tion of the collective qubit and phononic excitation even 
with a small number of ions. The features can be easily 
measured by local laser adressing. Compared to atoms 
in optical cavities, our implementation is easier to ma- 
nipulate, as all parameters can be tuned by changing the 
trap frequency, laser detuning and intensity. Addition- 
ally, the system can be extended by adding impurities of 
ions with different mass to the ion crystals, which allows 
for simpler addressing of the radial phonon modes and a 
separation of coexistent phases. 
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